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THERMODYNAMIC  MODEL  OF  CREEP  AT 
CONSTANT  STRESSES  AND  CONSTANT  STRAIN  RATES 

Anatoly  M.  Fish 


INTRODUCTION 

It  is  known  that  ductile  materials  under  applied 
constant  stress  experience  creep  deformation;  failure 
occurs  when  the  deformation  reaches  a  certain  value. 
Strain  rates  in  such  tests  are  functions  of  stress,  time, 
temperature,  and  other  factors.  If  the  tests  are  per¬ 
formed  so  that  the  strain  rate  remains  constant,  the 
stress  becomes  a  function  of  the  strain  rate  and  tem¬ 
perature.  It  is  obvious  that  the  physical-mechanical 
properties  of  the  material,  and  consequently  the 
rheological  parameters  of  the  creep  model  selected, 
cannot  and  must  not  depend  upon  the  loading  regime. 

Until  now  there  was  no  unified  creep  model  capable 
of  describing  a  material’s  behavior  under  both  loading 
regimes  in  the  wide  spectrum  of  stress,  time,  and  tem¬ 
perature.  The  strength  characteristics  of  viscous  mate¬ 
rials  (soils)  were  determined  from  constant  strain  rate 
(CSR)  tests,  while  their  deformation  characteristics 
were  determined  from  creep  (CS)  tests  at  constant 
stress.  However,  the  deformation  characteristics  ob¬ 
tained  from  CSR  tests  did  not  correspond  to  those 
from  CS  tests.  The  absence  of  adequate  information 
on  the  basic  stress-strain  relationship  during  creep 
affected  the  accuracy  and  reliability  of  solutions  to 
practical  problems. 

In  this  paper  for  the  first  time  the  entire  creep  proc¬ 
ess,  including  primary,  secondary,  and  tertiary  creep, 
and  failure,  for  both  constant  stress  tests  and  constant 
strain  rate  tests,  is  described  in  the  form  of  a  unified 
constitutive  equation  and  unified  failure  criteria.  Re¬ 
gardless  of  the  loading  regime  the  mechanical  proper¬ 
ties  of  the  material  are  described  by  the  same  param- 
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Figure  l.  The  creep  model. 

eters.  Hence,  information  obtained  from  one  type  of 
test  can  be  used  to  predict  a  material’s  behavior  under 
another  type  of  loading. 

This  creep  model  (Fig.  1 )  has  been  developed  based 
upon  several  theoretical  concepts: 

•  The  Eyring-Frenkcl  Theory  of  Rate  Processes 
(RPT) 

•  The  thermodynamic  nature  of  the  strength  of 
solids 

•  The  unity  of  the  process  of  deformation  and 
time-dependent  failure 


n’V.n’.V 


•  The  principle  of  the  linear  summation  of  de¬ 
fects  in  solids 

•  The  concept  of  normalized  strain,  strain  rate, 
stress,  and  time 

A  detailed  derivation  of  the  principal  relation¬ 
ships  of  the  model -that  is,  the  constitutive  equation, 
the  flow  equation,  and  the  failure  criteria  based  up¬ 
on  these  concepts-is  given  in  Fish  (1983b).  An 
approximate  derivation  of  the  model  is  given  in  Fish 
(1976a,  1980). 


The  creep  process  is  assumed  to  be  isothermal,  and 
volumetric  and  instantaneous  strains  are  assumed  to 
be  small. 

Failure  criteria 

When  t  =  tm,r=  1 ,  and /(I )  =  0,  the  strain  rate 
reaches  the  minimum  (em )  (points  m  in  Figs.  2  and  3), 
and  failure  occurs.  The  failure  state  is  characterized 
by  the  thermodynamic  failure  criterion  (Fish  1983b) 

AS  =  0  .  (4) 


PRINCIPAL  RELATIONSHIPS 
Constitutive  equation 

A  unified  constitutive  equation  describing  the  en¬ 
tire  creep  process  (including  primary,  secondary,  and 
tertiary  creep)  and  failure  at  constant  stress  has  been 
derived  by  Fish  (1976a;  1979;  1980;  1982a, b;  1983a, 
b): 

.  -kT  ,  E.  AS.  o.n+m  ... 
<-CFexp(-^).xpT<-)  (l| 

where  C,n>  0,  and  m  >  1  are  dimensionless 
parameters  independent  of  temperature 
aQ  =  temperature-dependent,  ultimate  strength 
of  soil,  MPa  (kg/cm2 ) 

E  -  activation  energy,  kJ/mole  (kcal/mole) 
k  =  Boltzmann’s  constant,  1.38x  10“23  J/°K 
h  =  Planck’s  constant,  6.63x  10'34  Js 
R  =  gas  constant,  8.31  x  10"3  kJ/mole  x  °K 
T  =  absolute  temperature,  °K 
kT/h  =  frequency  of  vibration  of  elementary  par¬ 
ticles  around  their  equilibrium  positions 
AS  =  change  of  entropy  (Fish  1983b) 


AS_AS(t)_ 


-Wt) 


(2) 


5  =  a  dimensionless  parameter  independent  of 
temperature 

fir)  =  a  function  of  normalized  (dimensionless) 

time 


/(t)  =  t  -  In  t  - 1 


(3) 


Figure  2.  Creep  curves  at  constant  stresses  ( CS)  (o  =  Const. ). 
Data  from  Zhu  and  Carbee  ( 1 983b).  m  are  the  inflection 
points,  ejj,  is  the  creep  strain  at  failure,  and  tm  is  the  time 
to  failure. 


and 


tm  =  the  time  to  failure. 


Figure  3.  Log  e  vs  log  t  for  constant  stress  f  CS)  tests.  Data 
from  Zhu  and  Carbee  ( 1 983a).  m  are  the  inflection  points, 
em  is  the  minimum  strain  rate,  tm  is  the  time  to  Jailure,  e, 
is  the  initial  strain  rate,  and  t,  is  the  initial  time. 
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Time  to  failure  tm  is  linked  to  applied  stress  a  by 
the  relationship  derived  by  Fish  (1976a,b;  1979;  1980; 
1982a,b;  1983a,b): 

n 

/m=/o(— )  for  o<a  (5) 


where  r0  is  the  mean  duration  of  the  “settled  life”  of 
an  elementary  particle  in  a  position  of  equilibrium 
(Frenkel  1947): 

h  E/RT 

l°*kfe  (6) 

and  hlkT  is  the  period  of  the  vibrations  of  the  elemen¬ 
tary  particles  around  their  equilibrium  positions. 

It  follows  from  eq  5  that  of  all  the  values  of  failure 
stress  amax  obtained  from  CSR  tests  (Fig.  4),  only 
one  stress,  o0,  causes  failure  at  time  ro.  The  basic 
postulate  of  this  author  is  that  time  tQ  is  equal  to 
Frenkel’s  relaxation  time  (eq  6).  The  stress  at  time 
rm  =  ta  is  defined  as  the  instantaneous  strength  of  the 
soil  a0. 

The  detailed  derivation  of  failure  criterion  eq  5 
and  the  physical  sense  of  the  stress  ratio  are  given  in 
Fish  (1983b),  where  it  was  shown  that 


m  /  o»  1 

r°  T"(v)  V“ 


where  r0  is  the  normalized  time  to  failure  and  w  and 
p  are  the  thermodynamic  and  mathematical  probabil¬ 
ities,  respectively,  of  the  concentration  of  m  defects 
(microcracks)  in  the  unit  volume  of  the  material 
(soil).  According  to  eq  7,  time  to  failure  tm  is  in  fact 
the  product 
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Figure  4.  Stress-strain  relationship  at  constant  strain 
rates  (CSR).  Data  from  Zhu  and  Carbee  (1983b).  m 
are  the  inflection  points,  em  is  th*  failure  strain,  and 
o  is  the  peak  vabe  of  stress. 


Since  stress  a  <  o0,  to  changes  from  unity  up  to  in¬ 
finity,  so  the  time  to  failure  cannot  be  smaller  than 
time  ta. 

If  the  magnitudes  of  tQ  obtained  from  tests  are 
smaller  than  those  calculated  by  eq  6  using  the  appro¬ 
priate  values  of  the  activation  energy,  the  loading  re¬ 
gime  is  considered  to  be  dynamic  (i.e.  an  impact  load¬ 
ing  regime)  rather  than  static.  At  the  same  time,  at 
high  strain  rates  (Sayles  and  Epanchin  1966,  Gold 
1973,  Bragg  and  Andersland  1980)  the  times  tQ  may 
vary  while  the  magnitude  of  oa  remains  unchanged. 
However,  the  times  t0  are  small  and,  for  practical  pur¬ 
poses,  such  variations  in  tQ  are  insignificant. 

Note  that  by  definition  the  logarithm  of  w  is,  with 
accuracy  up  to  a  constant,  the  entropy,  i.e. 

.  .  .  o0 

—  =  In  w  =  In  r0=m  In  (-).  (9) 

reflects  the  change  of  the  entropy  caused  by  the 
application  of  stress.  Equation  9  is  the  thermodynam¬ 
ic  failure  criterion  for  constant  stress  (CS)  tests.  For 
constant  strain  rate  (CSR)  tests,  when  o-  o  max,  the 
thermodynamic  failure  criterion  is 


4}  |  ,  .  , 

T  =  1"  =  m  In  (- )  . 

^  O  may 


When  o  =  omax  =  0o,  eqs  9  and  9a  yield  the  failure 
criterion 

A$  =  0. 

The  validity  of  eqs  9  and  9a  is  shown  below  (Figs.  5 
and  6)  using  test  data  on  creep  of  frozen  soil  at  con¬ 
stant  stresses  and  at  constant  strain  rates. 
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Figure  6.  Determination  of  parameters  m  and 
E(t„). 

Secondary  creep:  Flow  equations 

When  t  =  tm ,  t  =  1 ,  and  /(I )  =  0,  then  the  strain 
rate  e  becomes  em .  According  to  eq  1 ,  there  can  be 
three  types  of  relationships  between  the  minimum 
strain  rate  em  and  the  applied  stress  a  with  respect 
to  the  value  of  parameter  n. 

General  case  (0  <n¥=  1 ) 

In  this  case,  AS  =  0  and  eq  1  becomes  (Fish  1 979; 
1980;  1982a,b) 


-  kT  F  a  " 

tm=c~h  exP  (-jf) 

Combining  eqs  5, 6,  and  10, 


C  ,  a  "+m  C  ,o” 

em = r  =  r  (r> 

o  uo  ‘m  uo 


Denoting C  =  e*  Cand 


On  n+tn  n  - 


(10) 


(11) 


eq  1 1  takes  the  form 


•  _  C  -  Const- 

cm  “  t  t 

‘r  r 


(12) 


(13) 


where  the  time  to  rupture,  tT  >  tm,  is  the  time  to  the 
complete  destruction  of  the  sample  (not  shown  in 
Figs.  1, 2,  and  3).  The  relationship  in  eq  13  has  been 
confirmed  by  test  data  for  a  great  variety  of  materials 
(Regel  et  al.  1974). 

Note  that  stress  oa  in  eqs  1  through  1 2  is  not 


merely  a  reference  unit  to  balance  the  dimensions  of 
stress,  but  is  a  definite  strength  characteristic  depend' 
ent  upon  temperature  and  determined  from  constant 
strain  rate  (CSR)  tests  at  a  definite  (non-arbitrary) 
rate  of  deformation.  The  mathematical  structure  of 
the  stress  ratio  (o/a0)  is  identical  to  that  giver’  by 
Nadai  (1950),  Rabotnov  (1966),  Ladanyi  (1972), 
and  others.  From  the  physical  viewpoint,  parameter 
oQ  may  be  considered  to  be  the  strength  of  the  mate¬ 
rial’s  crystalline  lattice  (Garofalo  1965)  or  of  the  soil 
structure. 

Generalized  Newtonian  law 
of  viscous  flow  (n  =  l ) 

If  n  =  1 ,  then  denoting,  as  suggested  by  Fish  (1976a, 
1979, 1983b), 


E°  =  -^  and  rj  =  E0t0 
and  combining  eqs  1 1  and  14, 


•  _  o  ,  a . 1,1  _  a 

£m~EjnW  ~ET 


(14) 


(15) 


or 


.  a  ro  a  =o^  a 


(16) 


where  p  and  u>  are  the  mathematical  and  thermal  prob¬ 
abilities,  respectively.  Equation  1 6  can  be  considered 
as  a  generalized  law  of  Newtonian  viscous  flow.  When 
n  =  0  and  m  -  1  (see  below),  eq  16  transforms  into 
Newton’s  law  of  viscosity. 


.  o 

=  — 


(17) 


and  parameters  t0,E0,  and  rj  gain  a  physical  sense  of 
Frenkel’s  relaxation  time,  a  deformation  modulus,  and 
a  viscosity  coefficient,  respectively  (Frenkel  1947). 

Equations  14  and  1 5  establish  the  strongest  links 
between  the  deformation  modulus  of  soil,  the  viscos¬ 
ity  coefficient,  soil  strength  and,  as  will  be  shown,  the 
critical  failure  strain. 

Particular  case  (n  =  0) 

The  case  when  n  =  0  implies  that  the  slope  of  the 
plot  (log  em  vs  log  /m)  is  equal  to  -1  (Figs.  3,  7,  and 
8),i.e. 
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Figure  7.  Rheological  curve. 


various  materials  (Ladanyi  1972)  and  for  frozen  soil 
(Zhu  and  Carbee  1983a)  shows  that  the  deviation 
from  n  =  0  is  not  substantial.  From  tests  on  frozen 
Fairbanks  silt,  the  value  of  n  is  found  to  be  within  the 
limits  0  <  n  <  1  over  a  wide  spectrum  of  stresses  and 
temperatures.  Therefore,  when  comparing  creep  at 
constant  stress  and  constant  strain  rates,  the  value  of 
n  can  be  ignored. 


CREEP  AT  CONSTANT  STRESS 
(or  =  CONST.) 

Creep  model 

The  principal  concept  of  this  creep  model  is  that 
deformation  and  failure  compose  a  unified  process 
that  can  be  described  by  a  unified  constitutive  equa¬ 
tion  with  only  one  parameter,  6. 

Combining  eqs  1 , 2,  and  10,  we  obtain  (Fish  1983b, 
Fish  and  Assur  1983), 


(19) 


where,  as  a  first  approximation,  the  function 
/(t)  =  t  -  In  t  -  1  and  t  =  t/tm  . 


The  creep  process  can  therefore  be  characterized 
by  the  values  of  normalized  time  r,  strain  rate  |  =  i/em 
(see  eq  22  below),  and  the  corresponding  value  of 
function  /(r): 

Primary  creep  r  <  1  f(r)  >0  |  >  1 

Secondary  creep  r  =  1  /(l)  =0  \  =  1 

Tertiary  creep  t>  1  /(r)>0  |>1 


Figure  8.  Determination  of  parameters  E(t  0)and  (n+m). 

This  concept  is  supported  by  Ladanyi  (1972), 

Mitchell  (1976),  Assur  (1980),  and  partially  by  Mellor 
and  Cole  (1983),  and  others.  It  is  assumed  that  in 
short-term  creep  the  inflection  points  of  the  transition 
from  primary  creep  to  tertiary  creep  are  clearly  distin¬ 
guishable.  In  practice,  however,  it  can  be  difficult  to 
find  the  precise  locations  of  inflection  points  m  and 
to  determine  accurately  times  to  failure  tm .  For  short¬ 
term  creep,  slight  variations  of  the  absolute  values  of 
tm  (and  n)  do  not  affect  substantially  the  results  of 
the  final  creep  strain  calculations. 

It  is  also  possible  (Regel  et  at.  1974)  that  the  slope 
of  the  plot  (log  em  vs  log  tm )  is  not  equal  to  -1  (n  ¥=  0), 
as  described  by  eqs  1 0  and  1 3.  Analysis  of  test  data  for 


Equation  19  implies  that  the  duration  of  secondary 
creep  is  zero  (point  m  in  Fig.  1 ).  In  developing  this 
model  it  was  assumed  that  the  major  cause  of  deforma¬ 
tion  is  the  failure  process,  which  is  initiated  at  the  very 
start  of  deformation  as  soon  as  the  constant  load,  how¬ 
ever  small,  is  applied  to  the  sample. 

Equation  19  belongs  to  the  general  family  of  creep 
models  derived  from  various  assumptions  by  Zaretskii 
and  Vialov  (1971 )  and  Vialov  et  al.  (1973)  and  by 
Fish  (1976a,  1979, 1980).  The  simplest  conceptual 
approach  leading  to  eq  19  was  given  by  Assur  (1980). 
The  relationship  of  this  equation  to  primary  creep 
models  was  considered  by  Fish  and  Sayles  ( 1981 ). 

Generalized  model  of  viscous  flow 

If  n  =  I,  combining  eqs  15  and  19  (Fish  1979)gives 
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Figure  9.  Log  e  vs  normalized  time  t  for  constant 
stress  (CS)  tests,  m  are  inflection  points. 
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Figure  10.  Normalized  strain  rates  e/em  vs  normal 
ized  time  t  for  constant  stress  (CS)  tests,  m  is  the 
inflection  point. 
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Figure  11.  Normalized  stress  (omaxj of”  vs  nor- 
malized  time  t  for  constant  strain  rate  (CSR) 
tests,  m  is  the  inflection  point. 


(20) 


10, 


#  =  _°_ 

n(r)' 

where  the  function 

V(r)  =  E0tm  e-6f(r>  (21) 

is  the  coefficient  of  viscosity,  which  is  dependent  on 
stress,  temperature,  and  the  normalized  time  t. 

Normalized  creep  strain  rate 

According  to  eq  19,  the  behavior  of  various  mate¬ 
rials  during  creep,  regardless  of  stress,  temperature, 
or  time,  can  be  described  by  a  unified  constitutive 
equation: 


(22) 


where  |  is  the  normalized  (dimensionless)  strain  rate 
(Figs.  9  and  10).  Thus,  all  the  specifics  of  a  mate¬ 
rial’s  behavior  with  respect  to  time  are  expressed  by 
the  magnitude  of  parameter  5 : 


s‘”7w“ 


(23) 


Figure  12.  Determination  of  parameter  8. 

where  r,  =  r,/rm ,  t,  is  the  initial  (or  any  selected) 
time,  and  e,  is  the  corresponding  strain  rate  (Fig.  3). 

Creep  strain  ( o  =  Const.) 

General  case  (0  <  n  #  1) 

The  tabulated  integral  of  eq  19  for  8  =£  1  (Fish 
1983b,  Fish  and  Assur  1983)  gives  an  equation  for 
creep  strain  ec : 
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The  material’s  behavior,  with  respect  to  stress  and 
temperature,  is  characterized  by  the  absolute  values 
of  parameters  t0,o0,  m,  and  n,  which  compose  flow 
equation  eq  10  and  failure  criterion  eq  5. 

Note  that  from  the  physical  viewpoint  the  change 
in  the  normalized  strain  rate  over  time  reflects  the 
change  of  the  entropy  increment,  i.e. 


|(r)  =  exp  y  ■ 


(24) 


This  relationship  is  shown  in  Figures  1 0  and  1 1 . 


At  secondary  creep,  when  t  =  t  ,  r  =  1 ,  and 


/(I )  =  0,  then  regardless  of  the  material  type,  stress, 
or  temperature  eqs  22  and  24  yield  the  failure  crite¬ 
rion 


which  corresponds  to  eq  4. 

Parameter  8  can  be  found  from  Figure  1 2  as  a 
slope  of  the  straight  line,  or  it  may  be  calculated  at 
any  moment  of  time: 


6  = 


In 


t,  -  In  t,  -  1 


(25) 


e  t  \l> 


(26) 


or 


ec  =  e  t 

m  m 


^  t  e 


Sf(T) 


(27) 


where  is  the  integration  factor.  The  instantaneous 
strain  is  assumed  to  be  small  and  can  be  ignored. 
Combining  eqs  1 1  and  27, 


ec=C(— )  \jtre 


Sf(r) 


(28) 


and 


5  r 

v'~Xl  1+X  (1+X)  (2+X) 


(St)2 


-...]  (29) 


in  which  X  =  1  -  5. 

Experiments  have  shown  that  for  frozen  soils, 

0  <  6  <  1 .  In  general,  six  terms  of  the  series  (eq  29) 
provide  sufficiently  accurate  calculations  of  creep 
strain.  For  preliminary  analyses,  when  0  <  8  <  0.3 
and  t  <  rm ,  the  calculations  can  be  simplified  by  re¬ 
taining  the  first  term  of  the  series  (Fish  1982a,  1983a). 
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Sfra/n  at  failure 

When  t  =  tm, t  =  1, and/(l )  =  0,  then  eqs  26,  27, 
and  28  yield 


ec  =  e  f 

m  mm  v  m 


-  n  » 

£<--  =C(-)  0 

m  v  m 

vo 


where  is  creep  strain  at  failure  and  integration  con¬ 
stant  \pm  is  given  by 

=  X  11  "  (HX)  +  (1+X)  (2+xj  '  1  (32) 

where  X  =  1  -  6.  Since  a0  is  temperature-dependent,  it 
follows  that  in  the  general  case  the  magnitude  of  strain 
at  the  moment  of  failure  is  a  function  of  stress  and 
temperature. 

Particular  case  (n  =  l ) 

If  n  =  1 ,  we  can  combine  eqs  1 5  and  27  to  obtain 


ec  ^TeSf(r) 

ho 


When  t  =  tm ,  r  =  1 ,  and  /(I )  =  0,  eq  33  becomes 


ec  =  ~  Jj 

m  p  ym  • 


Equation  34  establishes  an  important  linear  relation¬ 
ship,  similar  to  that  of  Hooke’s  law,  between  the  ap¬ 
plied  stress  and  the  strain  at  failure.  In  practice,  how¬ 
ever,  such  a  dependency  may  not  always  be  recogniz¬ 
able  because  of  variations  in  the  soil  properties,  which 
are  expressed  in  variations  of  6,  Ea,  and  n. 

Particular  case  (n  =  0) 

If  n  =  0,  the  calculation  of  strain  by  eq  28  is  simpli¬ 
fied,  i.e. 

ec  =C4/  re{f(T)  (35) 

When  t  =  tm ,  t  =  I ,  and  /( 1 )  =  0,  eq  35  yields 

6m  =  &  tf'm  =  Const-  =  et  (36) 

where  ef  is  the  failure  strain. 

Equation  36  indicates  that  failure  strain  ef  is  con¬ 
stant  regardless  of  stress  and  temperature  and  can  be 


considered  as  the  failure  criterion  for  frozen  soil  and 
ice  (Ladanyi  1972,  Assur  1980,  Mellor  and  Cole  1982, 
Zhu  and  Carbee  1983a). 

A formalized  creep  curve 

Creep  curves  obtained  for  various  materials,  regard¬ 
less  of  the  magnitude  of  applied  stress  and  tempera¬ 
ture,  can  be  superposed  on  one  universal  plot,  as 
shown  in  Figure  1 ,  and  they  can  be  analyzed  and  de¬ 
scribed  by  a  single  equation.  Combining  eqs  27  and 
30, 

£  =  !L  =  ±  re6*W  (37) 

ec  r  m 


where  if  is  the  normalized  creep  strain.  At  secondary 
creep,  when  t  =  1  and  /( 1 )  =  0,  regardless  of  the  mate¬ 
rial  type,  stress,  or  temperature,  the  failure  criterion 
takes  the  form 

where  ifm  is  the  normalized  failure  strain. 

In  an  approximate  form  eq  37  yields 

$=.£.*- r  r)  (38) 


where  e£,  is  given  by  eqs  30, 3 1 , 34,  or  36,  in  which 
*m~l/(l-6)*.  *<1  (32a) 

is  obtained  by  retaining  the  first  two  terms  in  eq  32, 

/i(r)=r-61nr/6,  -1,  and  (3a) 

S^S+O-S)*-^-!-  -  ±  (32b) 

y  m  y  m 


Equation  38  combines  both  the  simplicity  and  the 
satisfactory  accuracy  of  the  creep  strain  calculations, 
and  retains  the  basic  assumptions,  the  parameters  of 
the  model,  and  their  physical  sense. 


CREEP  AT  CONSTANT  STRAIN 
RATE  (i  =  CONST.) 

Stress-strain  relationship 

If  the  tests  arc  performed  at  a  constant  strain  rate 
(e  =  =  Const.),  stress  o  in  eqs  1  and  9  becomes  a 

function  of  the  strain  rate  and  the  strain.  For  CSR 
tests,  assuming  that  strains  are  small,  i.e.. 


and 


e  =  et, 


The  rheological  curve  obtained  from  the  CSR  tests  is 
presented  in  Figure  7. 
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(39) 


Combining  eqs  10  and  39  and  taking  into  account 
that  when  e  =  em  (Figs.  1  and  4)  the  stress  reaches 
a  peak  (maximum): 
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At  the  same  time,  we  can  combine  failure  criterion 
eq  5  and  eq  39: 


Normalized  stress-strain  curve 

Thus,  three  principal  relationships  have  been  ob¬ 
tained  that  compose  the  creep  model :  the  flow  equa¬ 
tion  (eq  44a),  the  failure  criterion  (eq  44),  and  the 
stress-strain  relationship  (eq  43)  are  sufficient  to  de¬ 
scribe  the  entire  creep  process  during  constant  strain 
rate  tests. 

Determining  the  parameters  is  substantially  simpli¬ 
fied  if  all  the  test  data  are  presented  in  normalized 
form,  as  in  Figure  1 .  Thus,  denoting  a  =  o/a  max  gives 

d  =—  =e-sf(T)lm  (45) 

°max 


or 


t  °o(~) 


(41) 
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g  Sf(r)/m 


(46) 


Comparing  eqs  40  and  41 ,  we  see  that  the  simplest 
way  to  relate  CSR  tests  and  CS  tests  is  to  assume 
that  the  strains  are  small  and  that 


We  can  conclude  that  in  CSR  tests  deformation  and 
failure  also  compose  a  unified  process  that  can  be  de¬ 
scribed  by  a  single  equation  with  only  one  parameter, 
8/m: 


n-  0  and  C  =  em .  (42) 

Equation  41  is  in  agreement  with  the  failure  criteria 
for  CSR  tests  suggested  by  Sayles  and  Epanchin 
(1966),  Ladanyi  (1972),  Bragg  and  Andersland 
(1980),  Mellor  and  Cole  (1983),  and  others. 

Combining  eqs  10, 19, 41 ,  and  42,  a  constitutive 
equation  and  a  failure  criterion  for  tests  at  constant 
strain  rates  take  the  forms 


8  _  ln(pm J°) 

m  ~  At) 

where  /(t)  given  by  eq  3  is 

/(r)  =  7  -  In  r  -  1  and  t  =  t/tm 


(46a) 


Note  that  eq  45  can  be  rewritten  in  the  form 


o  =  o(t)  =  omax  e'bft'T)lm  (43 ) 
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respectively.  Time-dependent  failure  in  CSR  tests 
is  shown  in  Figure  5.  Since  in  CSR  tests  e  =  em , 
eq  44  can  be  presented  in  the  form  of  a  flow  equa¬ 
tion  similar  to  eq  1 8 : 


(44a) 


which  is  in  agreement  with  eqs  2, 9,  and  24.  This  re¬ 
lationship  with  respect  to  t  is  shown  in  Figure  1 1 . 

By  analogy  with  creep  at  constant  stress,  the  defor¬ 
mation  process  at  a  constant  strain  rate  can  be  charac¬ 
terized  by  the  values  of  the  normalized  time  (strain)  t, 
stress  a,  and  the  corresponding  value  of  function /(t): 


Primary  creep 
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/(7)>  0 
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Secondary  creep 

7  =  1 

/(1>  =  0 

0  =  1 

Tertiary  creep 

7  >  1 

/(7)>  0 
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As  in  the  case  of  short-term  creep  at  constant  stress, 
the  duration  of  secondary  creep  in  CSR  tests  is  also 
assumed  to  be  zero  (points  m  in  Figs.  1  and  4  and 
Fig.  13). 

When  time  /  =  /m,r  =  l,and/(l)  =  0,  then  re¬ 
gardless  of  the  material  type,  strain  rate,  or  tempera¬ 
ture  eq  45  yields  the  failure  criteria: 

a  =  1  and  e  =  1 

equivalent  to  eq  4.  The  normalized  stress-strain 
curve  constructed  according  to  eq  45  is  shown  in 
Figure  1 . 

Secant  modulus 

The  stress-strain  relationship  for  CSR  tests  can 
also  be  expressed  in  the  terms  of  a  secant  modulus 
(Fig.  14): 

a  =  eE(r)  (48) 

where 

E(e)  =  £(r)  =  e-«f(r)/m  (49) 

is  the  secant  modulus.  At  the  moment  of  failure, 

9  =  am«*'  *  *  em  •  and  *(»■)  ■  Em  -  '-e- 

°  max  -  em  ^  m  •  (^0) 

Combining  eqs  49  and  50,  one  obtains 
£ 

E( T)  =  —  e-6f(T)/m  .  (51) 
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Figure  13.  Determination  of  ratio  6/m. 


E,  Strain 


Figure  14.  Diagram  for  calculation  of  se¬ 
cant  modulus.  Em  is  the  secant  modulus 
at  failure. 

Note  that  eq  51  is  nonlinear,  similar  to  the  modu¬ 
lus  of  the  nonlinear  theory  of  elasticity.  The  magni¬ 
tude  of  the  secant  modulus  F(t)  approaches  infinity 
when  t  -*•  0  or  e  -*■  0.  For  solutions  of  boundary 
problems,  however,  such  a  deficiency  of  the  model 
is  insignificant. 

Stress/strain/strain  rate  at  failure 

The  relationship  between  stress,  strain,  and  the 
applied  strain  rate  at  failure  when  t  =  tm ,  r  =  1 ,  and 
e  =  e,,,  for  two  consecutive  tests  is,  according  to  eq 
41,  described  by 


i 


where  0maXl  and  em  are  the  peak  (maximum) 
stress  and  its  corresponding  strain,  respectively,  ob¬ 
tained  at  strain  rate  e  i ;  and  omaxj  and  emj  are  the 
peak  (maximum)  stress  and  its  corresponding  strain, 
respectively,  obtained  at  strain  rate  e2 . 

If  emi  =  =  Const.,  eq  52  is  simplified: 

CTmaxi  ~  ‘ 

Equation  53  is  in  agreement  with  a  known  empirical 
relationship  widely  used  in  describing  creep  in  metals 
(Nadai  1950,  Rabotnov  1966). 

TEST  DATA 

The  results  of  short-term,  unconfined  compression 


1 


creep  tests  performed  by  Zhu  and  Carbee  (1983a,b) 
on  Fairbanks  silt  (FS)  at  constant  stress  and  at  con¬ 
stant  strain  rates  are  presented  in  Figures  2, 3,  and  4 
and  Tables  1  and  2.  The  tests  were  conducted  at  -5°C 
using  remolded  soil  specimens  with  a  dry  density  of 
from  1 .16  to  1 .24  g/cm3  and  a  water  content  from 
40%  to  44.8%,  respectively.  The  specimens  were 
vacuum-saturated  and  quick-frozen  in  an  open  system. 
The  CS  creep  tests  were  carried  out  using  a  test  appa¬ 
ratus  developed  by  Sayles  (1968).  The  CSR  tests 
were  performed  on  a  universal  loading  machine. 

Preliminary  analysis 

The  unity  of  the  deformation  and  failure  processes 
under  two  different  loading  regimes  -  constant  stress 
and  constant  strain  rate  —  will  now  be  demonstrated. 

Time-dependent  failure 

For  the  CS  creep  tests,  the  dependency  shown  in 
Figure  5  is  based  upon  the  test  data  presented  in  Fig¬ 
ure  2  and  Table  1 .  The  times  to  failure  /  for  the 

m 

corresponding  stresses  a  were  determined  at  points  m 
where  the  strain  rates  are  the  minima  em .  In  the 
normalized  time  to  failure  r0,  tn  =  2.7x  10~2  s  (see 
Constant  stress  tests  below).  The  stress  ratios 
(o/o0)  were  calculated  assuming  that  the  highest 
stress  in  the  CSR  test  is  the  instantaneous  strength 
0O  =  121.2  kg/cm2  (Table  2). 


For  CSR  tests  (Fig.  4),  the  stress  ratios  (omax/a0) 
were  calculated  using  the  same  stress  o0  =  1 2 1 .2  kg/cm5 . 
The  times  to  failure  /m  presented  in  Table  2  were  de¬ 
termined  by  eq  39,  using  the  strains  em  and  strain  rates 
e  that  correspond  to  the  peak  stresses  (amax );  param¬ 
eter  t0  =  2.7x  10'2  s  was  calculated  by  eq  6  (see  Con¬ 
stant  strain  rate  tests  below). 

In  Figure  5  the  experimental  points  for  both  load¬ 
ing  regimes  are  superposed  on  one  curve,  indicating 
that  in  both  cases  time-dependent  failure  can  be  de¬ 
scribed  by  a  unified  failure  criterion.  A  similar  rela¬ 
tionship  between  the  two  loading  regimes  for  poly- 
crystalline  ice  was  shown  by  Mellor  and  Cole  (1983). 

Rheological  curve 

In  the  classical  sense,  the  rheological  curve  is  the 
relationship  between  the  secondary  creep  strain  rate 
and  applied  stress,  presented  on  logarithmic  coordi¬ 
nates.  It  is  more  convenient,  however,  to  present  the 
test  data  according  to  eqs  10  and  44a,  as  a  dependency 
between  the  minimum  strain  rate  emand  the  stress 
ratio  (o/o0)  for  CS  tests,  and  as  a  dependency  be¬ 
tween  the  applied  strain  rate  (e  =  em)  and  the  stress 
ratio  (omax/o0)  for  CSR  tests.  The  same  instantane¬ 
ous  strength  (oD  =  1 21 .2  kg/cm2 )  was  used  to  calcu¬ 
late  both  stress  ratios. 

The  rheological  curve  for  both  loading  regimes  is 
presented  in  Figure  7.  The  experimental  points  for 


Table  1 .  Constant  stress  tests  (o  =  Const.).  Data  from  Zhu  and 
Carbee  (1983a). 
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Table  2.  Constant  strain  rate  tests  (e  =  Const.).  Data  from  Zhu 
and  Carbee  (1983b). 
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both  loading  regimes  are  superposed  on  one  curve,  in¬ 
dicating  that  in  both  cases  the  flow  of  soil  at  second¬ 
ary  creep  can  be  described  by  a  unified  flow  equation. 
A  similar  conclusion  was  reached  by  Mellor  and  Cole 
(1983)  based  on  test  data  on  the  creep  of  polycrystal¬ 
line  ice. 

Normalized  strain  rates  and  stress 
ratios  vs  normalized  time 

Fish  (1976, 1979,  1980)  introduced  this  creep 
model,  including  the  concepts  of  normalized  time 
( t/tm ),  normalized  stress  ratios  (a/o0),  and  conse¬ 
quently  normalized  strain  rates  (e/em ),  because  the 
conventional  plot  (log  e  vs  log  t)  shown  in  Figure  3 
and  the  conventional  rheological  curve  (log  em  vs 
log  a),  which  are  usually  used  to  derive  creep  models, 
were  inadequate. 

Indeed,  we  cannot  take  the  logarithm  of  a  quantity 
that  has  dimensions.  When  we  do  deliberately,  to 
emphasize  the  differences  between  particular  tests, 
we  always  assume  that  the  values  of  the  strain  rates 
and  the  time  have  denominators  (referenced  param¬ 
eters)  that  render  the  strain  rates  and  the  times  dimen¬ 
sionless.  In  the  case  of  creep,  however,  these  refer¬ 
enced  parameters  become  functions  of  stress  and  tem¬ 
perature,  and  instead  of  a  single  plot  (log  e  vs  log  /) 
on  which  all  the  test  results  would  be  superposed  re¬ 
gardless  of  stress  and  temperature,  we  obtain  a  sepa¬ 
rate  family  of  curves  for  each  temperature. 

It  has  always  been  known  that  it  is  virtually  im¬ 
possible  to  derive  empirically  a  constitutive  equation 
for  a  creep  model  or  to  select  an  empirical  formula  that 
corresponds  to  the  test  data  with  parameters  that  are 
independent  of  stress  and  temperature.  Every  creep 
curve,  even  for  a  given  stress  and  temperature,  has  at 
least  a  slightly  different  shape  due  to  variations  in  the 
properties  of  the  soil  specimen.  The  use  of  normalized 
concepts  together  with  the  Rate  Process  Theory  makes 
it  possible  to  smooth  out  differences  in  the  specific  be¬ 
haviors  of  specimens  and  to  derive  constitutive  equa¬ 
tions  with  temperature-independent  parameters. 

The  normalized  strain  rates,  stress  ratios,  and  times 
presented  in  this  paper  have  a  definite  physical  sense: 
they  reflect  the  change  of  the  entropy  of  the  creep 
process. 

It  is  precisely  the  entropy  increments  AS  and  AS 
that  unify  the  two  loading  regimes.  To  prove  the 
validity  of  this  statement,  the  test  data  in  Figures  2 
and  3  must  be  superposed  on  the  test  data  in  Figure  4. 

The  relationship  (log  e  vs  log  t)  obtained  in  the  CS 
tests  presented  in  Figure  3  can  be  replotted  on  semi- 
logarithmic  coordinates  (log  e  vs  r),  as  shown  in  Fig¬ 
ure  9.  Then  the  ordinates  of  the  curves  are  expressed 
in  terms  of  a  normalized  strain  rate  (£  =  e/em ),  divid¬ 
ing  the  strain  rates  i  by  the  corresponding  minimum 
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strain  rates  em.  This  plot  is  shown  in  Figure  10.  We 
can  see  that  the  experimental  points  form  a  curve 
that  can  be  described  by  a  single  constitutive  equa¬ 
tion.  According  to  eqs  2  and  24,  function  |(t)  indi¬ 
cates  the  change  of  entropy  increment  A5(r)  with 
normalized  time. 

The  relationship  of  the  stress  ratio  (omax/o)m  vs 
the  normalized  time  r  for  m  =  7,  based  upon  data  ob¬ 
tained  from  CSR  tests,  is  shown  in  Figure  1 1 .  In  this 
plot  the  experimental  points  also  form  a  curve  that 
can  be  described  by  a  single  constitutive  equation. 
According  to  eq  47,  the  change  of  the  function 
(l/a)m  =  (amax  lo)m  with  time  (r)  reflects  the  change 
of  the  entropy  increment  AS(t)  with  normalized  time. 

Figures  1 0  and  1 1  are  drawn  to  the  same  scale, 
so  that  the  plots  for  the  two  loading  regimes  can  be 
superposed.  In  Figure  1 1  the  CS  test  data  at  o  =  20 
kg/cm2  is  superposed  on  the  CSR  test  data.  Note 
that  in  both  cases  the  curves  have  minima  when 
time  t  =  1 ,  confirming  the  validity  of  the  failure  cri¬ 
terion  eq  4. 

Constant  stress  tests  (o  =  Const.) 

Parameters  m  and  tJE) 

To  determine  parameters  m  and  tQ  according  to 
eq  5,  the  test  data  on  time-dependent  failure  in  Table 
1  are  presented  on  the  coordinates  shown  in  Figure  6: 

tmkT  £  o 

T\n(~—  )=  ^mrin(^).  (5a) 

The  slope  of  the  straight  line  gives  m^l .  The  inter¬ 
sect  with  the  ordinate  is  equal  to  the  ratio  E\R 
6.9x  103  (°K),  which  gives  the  value  of  the  activation 
energy  E  **  13.8  kcal/mole.  Time  to  failure  fmin  this 
calculation  is  in  seconds  and  T=  268.16°K.  Parameter 
ta  =  2.7x  10'2  s,  calculated  by  eq  6,  is  smaller  than 
time  tmfor  oD obtained  in  the  CSR  tests  (Table  2). 

The  absolute  values  of  tQ  and  m  are  closely  inter¬ 
connected.  In  the  case  of  unfrozen  soils,  or  if  the  tem¬ 
perature  of  the  frozen  soil  is  constant,  ta  and  m  may 
be  determined  directly  by  plotting  the  test  data  with 
eq  5b: 


In  tm  =  In  ro  +  m  In  (— ) 


or  they  may  be  calculated  with  eq  5c: 


In  (t—)  =  m  In  (~) 
o 


in  which  tn  corresponds  to  time  to  failure  tm  for 
stress  o„  (Table  2). 


Note  that  ta  and  m  depend  substantially  upon  tem¬ 
perature.  For  instance,  calculated  by  eq  Sb,  m  ranges 
from  approximately  3  to  13  if  the  temperature  of  the 
frozen  soil  changes  from  -1°C  to  -10°C.  The  magni¬ 
tudes  of  t0  also  vary  within  wide  limits.  Determina¬ 
tion  of  t0  and  m  by  means  of  eq  Sa  gives  more  stable 
values  for  these  parameters. 

Parameters  8  and  ipm 

According  to  eq  23,  the  magnitude  of  5  is  deter¬ 
mined  by  plotting  the  normalized  strain  rate  (e/em) 
vs  the  function /(r)  on  semilogarithmic  coordinates. 
This  plot  is  shown  in  Figure  1 2  for  the  creep  curve  at 
a  -  20  kg/cm2  for  which  the  slope  of  8  =  0 39.  For 
comparison  with  the  CSR  tests,  the  data  from  only 
one  CS  creep  curve  and  one  CSR  curve  were  plotted 
in  Figure  12.  It  can  be  shown,  however,  that  the  en¬ 
tire  family  of  creep  curves  plotted  in  this  manner  also 
forms  a  straight  line.  In  this  case,  however,  the  slope 
8  may  be  slightly  different. 

The  integration  coefficient  \pm  =  1 .295  was  calcu¬ 
lated  by  eq  32.  \jjm  and  8  can  also  be  determined  from 
eqs  30  and  32a  using  test  data  (Table  1 ). 

Parameters  Cand  (n+m) 

The  values  of  parameter  C  are  calculated  by  eq  36 
using  the  strains  ecm  given  in  Table  1 .  Thus,  for  the 
test  at  o  =  20  kg/cm2 ,  ej^  =  9.2%,  =  1 .295,  and 

C=  0.071. 

The  sum  of  (n+m)  is  determined  according  to  eq 
10  by  plotting  the  test  data  in  Table  1  on  the  coordi¬ 
nates,  as  shown  in  Figure  8: 

7’ln(rf)=f +  ("+"l)7'ln(~).  (10a) 

m 

The  intersect  with  the  ordinate  is  equal  to  the  ratio 
E/R  as  6.9x  103  (°K).  The  slope  of  the  straight  line 
is  equal  to  (n+m)  a=  6.8,  which  approximately  equals 
the  magnitude  of  m  a*  7  obtained  above.  Thus,  since 
by  definition  parameter  n  >  0  (eq  1 ),  we  can  conclude 
that  for  these  tests  n  =  0. 

For  unfrozen  soils,  or  if  the  temperature  of  the 
frozen  soil  is  constant,  C  and  the  sum  (n+m)  may  be 
determined  directly  by  plotting  the  test  data  with  eq 
11a: 

Mem/0)  =  lnC+(»+m)ln(^-)  (11a) 

where  t0  corresponds  to  time  to  failure  tm  of  the 
instantaneous  strength  aQ.  As  in  the  case  of  time- 
dependent  failure,  however,  t0  and  (n+m)  can  be  very 
sensitive  to  temperature  change. 

In  previous  publications  (Fish  1982a,b;  1983a), 


other  methods  for  determining  parameters  8,  m,  n, 
and  t0(E)  have  been  discussed. 

Constant  strain  rate  tests  (e  =  Const.) 

Parameters  m  and  tJE) 

Since  the  equation  for  time-dependent  failure  at 
constant  strain  rates  (eq  44)  is  identical  to  that  for 
creep  tests  (eq  5),  the  method  for  determining  param¬ 
eters  m,  tQ,  and  E  is  analogous  to  that  for  creep  tests. 
Thus,  in  Figure  6  the  test  data  presented  in  Table  2 
are  plotted  in  accordance  with  eq  44: 

/ emkT\  p  a 

rln\ir)’il*mrl"(^>-  <«b> 

In  Figure  6,  we  see  that  the  CSR  test  data  are 
superposed  on  the  same  straight  line  as  the  CS  creep 
test  data.  Therefore,  the  values  of  parameters  m,  tD, 
and  E  will  also  be  the  same. 

For  unfrozen  and  in  some  cases  frozen  soils,  as 
shown  above,  the  determination  of  m  and  tQ  may  be 
simplified,  using  eq  1  la  in  which  o  =  a max,  n  =  0, 
and  C  =  em . 

Parameters  8  and  C 

According  to  eq  47,  when  m  is  known,  the  CSR 
test  curves  can  be  superposed  on  the  CS  creep  curves. 
This  is  shown  in  Figure  1 2  for  two  tests:  at  a  =  20 
kg/cm2  and  at  i  =  1 .1  x  10~s  s'1 .  The  value  of  8  re¬ 
mains  unchanged.  From  eq  46  it  follows  that  pa¬ 
rameter  8  can  also  be  obtained  by  plotting  the  stress 
ratios  (o  max/o)  versus  the  corresponding  values  of 
function  /(r)  on  semilogarithmic  coordinates,  as  in 
Figure  13.  In  this  case  the  slope  of  the  straight  line 
is  equal  to  8/m  =  5.5 x  10~2 ,  which  gives  5  =  0.385 
form  =  7. 

Parameter  C  in  the  CSR  test  is  by  definition  (eq 
42)  equal  to  the  failure  strain  em  given  in  Table  2. 

THE  PRINCIPLE  OF  SUPERPOSITION 

We  can  now  formulate  the  Principle  of  Superposi¬ 
tion  of  CS  and  CSR  tests: 

Creep  curves  at  constant  stresses  and  stress- 
strain  curves  at  constant  strain  rates  can  be 
superposed  if  the  test  data  are  presented  as 
functions  of  the  change  of  entropy  over  nor¬ 
malized  time. 

This  principle  can  be  expressed  in  terms  of  normal 
ized  strain  rates  and  normalized  stress  ratios,  combin¬ 
ing  eqs  22,  24,  and  47 : 


(54) 


max 

a 


m 

) 


where  the  left  and  right  sides  of  eq  54  refer  to  the 
data  obtained  from  CSand  CSR  tests,  respectively. 

The  Principle  of  Superposition  can  be  extended  to 
the  area  of  creep  deformation  and  the  stress-strain 
relationship.  Thus,  in  Figure  15  the  creep  curves 
from  Figure  2  are  superposed  on  the  curve  calculated 
by  eq  37  using  5  =  0.39  and  i/*m=  1 .295.  In  Figure 
16  the  stress-strain  curves  from  Figure  4  are  super¬ 
posed  on  the  curve  calculated  by  eq  45  using  5  =  0.39 
and  m  =  7. 

The  creep  model  can  be  considered  valid  if  the 
families  of  both  the  creep  and  the  stress-strain  curves 
(Figs.  2  and  4)  can  be  superposed  at  all  stages  of 


creep  —  primary,  secondary,  and  tertiary.  In  Figure 
17,  the  ec  and  cutest  data  from  Figure  15  and  the 
(omax  /o)1"  ,est  data  from  Figure  16  (m  =  7)  are  re¬ 
plotted  on  the  coordinates: 

Creep  at  constant  stress  (o  =  Const.) 

ec  \1/ 

Y  = - ™  =T.  (37a) 

ec  \1/ 
m  Y 

Creep  at  constant  strain  rate  (e  =  Const.) 

,"  =  (0^,mTe^f(r)=T  (473) 

o 

In  both  cases  the  test  data  form  straight  lines  at 
angles  of  45°  to  the  axis  of  the  normalized  time. 


Figure  15.  Normalized  creep  curve. 
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Figure  16.  Normalized  stress-strain  curve. 
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Note  that  both  entropy  increments  AS  and  AS 
are  functions  of  stress  and  temperature.  They  are 
closely  interconnected,  and  their  relationship  pre¬ 
determines  the  direction  of  the  creep  process: 

•  If  stress  a  =  Const.,  then  AS  =  Const,  and,  since 
AS  =  /(t),  the  strain  rate  is  also  a  function  of  time 
(eq  19). 

•  If  strain  rate  i  =  Const,  and  AS  =  /(r),  then  the 
entropy  increment  and  the  stress  in  eq  9a  both  be¬ 
come  functions  of  time: 


t  ,57> 

as  described  by  eqs  43  and  47. 

•  If  AS  =  0  (secondary  creep)  and  the  minimum 
strain  rates  in  CS  tests  are  equal  to  the  strain  rates 
in  CSR  tests  (em  =  e),  then  the  entropy  increments 
AS  in  both  CS  and  CSR  tests  are  equal  as  well,  and 
so  are  the  stresses  a  =  o  .  The  reverse  is  also  true. 


Figure  /  7.  Superposition  of  creep  (CS)  and  con¬ 
stant  strain  rate  (CSR)  tests. 

Combining  eqs  37a  and  47a,  we  can  write  the  Princi¬ 
ple  of  Superposition  for  creep  strain  in  analytical 
form: 


n,  =  "max"1 


The  CS  and  CSR  test  data  are  in  fact  superposed 
in  Figure  17,  confirming  the  validity  of  the  mod¬ 
el.  The  slight  diversion  at  the  ends  of  the  tertiary 
creep  stages  can  be  attributed  to  changes  in  the 
properties  of  the  soil  specimens  and  to  the  change 
of  the  failure  mechanism  at  this  stage. 


THERMODYNAMIC  EQUATION  OF  CREEP 

It  has  been  shown  that  deformation  and  failure 
in  CS  and  CSR  tests  are,  in  fact,  a  unified  thermo- 
activated  process  in  which  the  dominant  role  belongs 
to  the  change  of  entropy.  Thus,  combining  eqs  1 , 2, 
9,  and  9a,  and  assuming  that  n  =  0,  we  obtain  a  ther¬ 
modynamic  equation  of  creep  for  both  CS  and  CSR 
loading  regimes: 


*  If  strain  rate  e  =  eD  =  Const,  and  AS  =/(t),  then 
stress  is  also  a  function  of  time,  o  =  /(r).  At  second¬ 
ary  creep  when  AS  =  0,  the  stress  becomes  omax  =  oD. 
(It  is  assumed  that  this  case  is  not  applicable  for  CS 
tests  where  a<  oa .) 

•  If  AS  =  AS,  the  strain  rate  becomes  (Fish  1976a, 
1980) 


.  -kT  .  E  ,  C 

eo  =  CTexp(-—  )  =  -, 


and  when  AS  =  AS  =  0,  the  corresponding  peak  (max 
imum)  stress  is,  as  defined,  die  instantaneous  strength 

V 

The  creep  process  for  both  CS  and  CSR  tests  can 
be  expressed  in  terms  of  the  free  energy  E  of  the 
thermal  activation.  Thus,  we  can  substitute  Boltz¬ 
mann’s  constant  k  in  eq  56  for  the  gas  constant 
R  =  kN,  where  N  is  Avogadro’s  number,  and  AS  = 

A SN,  AS  =  ASN,  and 

E  =  E-AST+  AST  .  (59) 

If  AS  =  0  and  o  =  omax  (secondary  creep), 

E  =  E  +  AST  =  E  +  RT\n  r„  .  (60) 
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When  AS  =  0  and  a  =  omax  =  oa  (primary  and  terti¬ 
ary  creep), 

E  =  E-&ST=E-RT8f(T).  (61) 

Combining  eqs  56  and  59,  the  thermodynamic 
equation  of  creep  transforms  into  the  Eyring-Frenkel 
(Rate  Process  Theory)  kinetic  equation: 


e  =  C—  e'EIHT 
C  h 


in  which  the  free  energy  A' and  consequently  the  strain 
rate  i  can  be  calculated  using  only  three  parameters, 
C,  A,  and  5.  The  methods  for  determining  these 
parameters  from  CS  and  CSR  tests  have  been  dis¬ 
cussed  above. 

Since  A S  and  AS  are  temperature-dependent  and 
AS  (and  AS  in  CSR  tests  as  well)  is  a  function  of 
time,  the  apparent  activation  energy  E  becomes  a 
function  of  stress,  time,  and  temperature  (Tietz  and 
Dorn  1956;  Zhurkov  and  Sanfirova  1958;  Andersland 
and  Akili  1967;  Mitchell  et  al.  1968;  Andersland  and 
AlNouri  1970;  Regel  et  al.  1974;  Vialov  1978:  and 
others). 


CONCLUSIONS 

A  unified  creep  model  has  been  developed,  which 
for  the  first  time  describes  the  entire  creep  process  - 
including  primary,  secondary,  and  tertiary  creep  - 
and  failure  for  both  constant  stress  (CS)  tests  and 
constant  strain  rate  (CSR)  tests. 

•  Deformation  and  failure  are  considered  as  a 
single  thermoactivated  process  and  are  de¬ 
scribed  by  a  unified  thermodynamic  equation 
of  state.  The  activation  energy  was  found  to 
be  a  material  constant  independent  of  the 
loading  regime. 

•  The  dominant  role  in  the  processes  of  defor¬ 
mation  and  failure  belongs  to  the  change  of 
entropy:  failure  occurs  when  the  entropy 
change  is  zero. 

•  The  stress-strain  state  of  the  material  with  no 
changes  of  entropy  is  characterized  by  the 
minimum  values  of  the  strain  rates  in  CS  tests 
and  the  maximum  (peak)  values  of  stress  in 
CSR  tests. 

The  Principle  of  Superposition  has  been  formu¬ 
lated,  according  to  which  creep  (e  vs  /)  curves  and 
stress-strain  (o  vs  e)  curves  can  be  superposed  if  the 
data  from  CS  and  CSR  tests  are  presented  as  func¬ 
tions  of  the  change  of  entropy  over  normalized  time. 


A  method  has  been  developed  to  determine  the 
parameters  of  the  model  so  that  creep  deformation 
and  the  stress-strain  relationship  of  a  ductile  mate¬ 
rial  such  as  soil  can  be  predicted  using  information 
from  either  type  of  test. 
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